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Abstract 



We study nonequilibrium steady states of the driven lattice gas with two particles, 
using the most general stochastic transition rules that satisfy the local detailed 
balance condition. We observe that i) the universal l/r d long range correlation may 
be found already in the two-particle models, but ii) the magnitude (or even the 
existence/absence) of the long range correlation depends crucially on the rule for 
transition rates. The latter is in stark contrast with equilibrium dynamics, where all 
rules give essentially the same results provided that the detailed balance condition 
is satisfied. 

Stochastic processes with discrete state space are often studied as convenient idealized 
models of various physical systems. As for time evolution at or close to equilibrium, it 
has been established [H E] that transition rates must satisfy the detailed balance condition 
in order for the system to obey a macroscopic symmetry known as reciprocity. But no 
such criteria are known for systems far from equilibrium. A standard convention (see, for 
example, [HI EE]) is to take transition rates satisfying the local detailed balance condition, 
which is a direct generalization of the detailed balance condition. It seems that an implicit 
assumption has been that the specific choice of transition rates does not affect physics in 
a serious manner. We here show that this is far from the case in driven nonequilibrium 
systems . 

We here study the driven lattice gas [3J EJ E], which is one of the standard models of 
nonequilibrium systems driven by an external force. We use the most general transition 
rates satisfying the local detailed balance condition, and obtain the exact steady states 
of the models when there are only two particles. We first drive a simple condition under 
which the steady state is the same as the equilibrium state. The condition is satisfied in a 
large class of rules, but not in some of the standard ones including the Metropolis and the 
heat bath rules. Then we turn to transition rates which do not satisfy the condition, and 
investigate nonequilibrium corrections to the steady states. We find that, in the dimensions 
d > 2, the two-particle models may exhibit the universal l/r d long range correlation, 
which is often regarded as an essential feature of nonequilibrium steady states. Moreover 
the magnitude of the long range correlation depends crucially on the choice of transition 
rates: it is extremely large in the Metropolis rule while it is absent in the exponential rule. 
Although there were some remarks [6J about rule dependence in the driven lattice gas, 
they were mainly quantitative and did not show the sharp implication on the long range 
correlation. As for the one- dimensional lattice gas driven by boundary conditions, it was 
shown recently [7j that the model with the exponential rule has a long-range correlation 
while that with the zero-range rule has no correlation. 
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Model: Let A = {-(L - l)/2, ...,(L— l)/2} d C Z d be the d-dimensional Lx ■■■ x L 
hypercubic lattice where d > 2. We impose periodic boundary conditions. There are 
two identical particles on A, whose positions are denoted by x G A and y G A. The 
configuration of the model is given by a pair (x, G A x A, where we do not exclude 
the possibility of x — y. (See the "standard DLG" section for the treatment of hard core 
repulsion.) The two particles interact with each other via a potential V(x — y) = V(y — x), 
and are driven by a constant (electric) field E = (E, 0, . . . , 0). 

We denote by U the set of 2d unit vectors. Let us define a stochastic dynamics by 
specifying transition rates for any (x, y) and 8 G U as 

c[(x, y)^(x + 8, y)] = $(V(x - y), V(x + 8 - y); E • 8), 

c[{x, y) - (x, y + 8)] = <S>(V{x - y), V(x - y - 8); E • 8), (1) 

where §(v,v';E) is a function (of three variables) satisfying the local detailed balance 
condition 

<$>{v,v';E) = e v - v ' +E $(v',v;-E). (2) 

This reduces to the standard detailed balance condition if E = 0. In many cases jSj one 
uses $ which can be written as $(t> , v'; E) = <f)(—v+v' — E) with a function <p{h) (of a single 
variable) satisfying <p{h) = e~ h <fi(— h). Among the standard choices are i) the exponential 
rule with (p(h) = e~ h ^ 2 , ii) the heat bath (or Kawasaki) rule with (p(h) = 2/(1 + e h ), and 
iii) the Metropolis rule with <p(h) — 1 if h < and <p(h) = e~ h if h > 0. 

We are interested in the properties of the steady state distribution (stationary measure) 
p X:V , which is the unique solution of 

J2{-Px, y c[(x,y) -> (x + 5,y)]-p Xt yc[(x,y) -> (x,y + 5)] 

c[(x + 8, y) (x, y)] + p x>y+s c[(x, y + 8) -> (x, y)] } = 0, (3) 

for any (x, y) G A x A. By using the translation invariance, one can write the steady 
state distribution as p X:V = e~ v<yX ~ yS} p x - y . Note that we have extracted the equilibrium 
distribution. We find from © that p z is determined by 

J2 [-Pz ~c(z^ z + 8)+ p z+s c(z + 8^ z)\ = 0, (4) 

for any z G A, where 

c{z -> z + 8) = e- y(z) ®( V ( Z ), V (z + 5 )> ( E • 5)°), ( 5 ) 

(T=±l 

are the effective transition rates. 

Exponential condition: If the effective transition rates satisfy 

c(z -> z + 8) = c(z + 8 -> 2), (6) 

for any z G A and 8 G W, then the solution of (@J) is simply given by p z = const. This 
means that the stationary distribution p x ^ y is the same as the equilibrium distribution. 



2 



Let us examine the condition ©. Let e\ = (1, 0, . . . , 0), and 14' be the set obtained 
by removing ±ei from U. For 8 G 14' (where we have E • 8 = 0) the condition (JBJ) is 
automatically satisfied because of the (local) detailed balance condition (J2J with = 0. 
For 8 = ±ei, the condition © reads e _,u {$(u, v'; E) + $(?;,?/; = e~ v '{$(v',v;E) + 

$(f',t>; — E)}, which, using ©, can be written as 

$(v,v';E) = e E $(v,v';-E), (7) 

for any v, v', and E. This exponential condition is satisfied if $(t>, t>'; .E) = e E ^ 2 ^(v,v') 
with any equilibrium transition rate \fr(v,v') satisfying the detailed balance ty(v,v') = 
e v-v \fr(v>^ v y Among the three standard rules we mentioned, only the exponential rule 
satisfies the exponential condition (J7J). 



Long-range correlation: We now turn to the rules where the exponential condition (J7|) 
is not satisfied, and investigate nonequilibrium corrections. For simplicity we only consider 
the case with on-site interaction described by the potential such that V(o) = v ^ and 
V(z) = if z 7^ o, where o = (0, . . . , 0) G A is the origin. The behavior of nonequilibrium 
corrections is essentially the same for more complicated interactions. 

Since v (z) = for z ^ o, one has c(z — > z ± e{) = s for any z G A unless z or z ± e\ 
is o. Similarly for any 8 G 14', one has c(z — » z + 8) = t unless z or z + 5 is o. The only 
transition rates which are not equal to s or t are c(o — > ±ei) = s', c(±ei — > o) = s", 
and c(o — > 8) — c(8 — > o) = t' where 8 G 14'. We assume s' 7^ s" since the exponential 
condition ((7j) is satisfied if s' = s" . 

Let us write the condition (jlj) for stationarity as J2 z 'eA{T z ,z' + D ZjZ i)p z i = for any 
2; G A. Here we defined the hopping matrix iT g)Z i) z>z i e \ by T Z)Z = —2s — 2(d — l)t, 
T Z)Z ± ei = s, T z 2+( 5 = t for all z G A and 8 G W, and T Z)Z / = otherwise. We also defined 
(D z> z')z, z >eA by D 0j0 = 2(s - s') + 2(d - l)(t - t'), D ei , ei = D_ ei ,_ ei = s - s", D M = * - t', 
-D ,± ei = s" — s, D± eit0 — s' — s, D 0y s = -D<5, — t' — t for all 5 G W, and -D^y = otherwise. 

We write the (unnormalized) steady state distribution as p z = 1 + ip z and study the 
behavior of the nonequilibrium correction ip z . Since Y^, z 'T ZiZ i — 0, the equation that 
determines ip z is 

^2 (T x>z i + D z>z ,) tp z , = -q z , (8) 
z'eA 

for any z G A, where the "charge" is defined by q z = Y^z'eA^z,z'- Here it is found that 
q Q = 2Q, q± ei = —Q, and q z = otherwise, where Q = s" — s' 7^ 0. 

Note that one has ipz = if v = or E = 0. Therefore, in the lowest order of 
perturbation in v and E, the nonequilibrium correction if) z is determined by the standard 
lattice Poisson equation J2 Z ' T Z)Z ' 4>z' = —Qz- Since the charge q z describes a quadrupole 
with the magnitude Q, the asymptotic long-distance behavior of ip z for d > 2 is readily 
obtained as 

r) 2 v 2 d 2 

(d-i)( gl y a )-E? =2 (^ 2 A) 

c ^ r/_. 9/_\ , 9 mi ' v y / 



{(^i 2 A) + E? =2 (^ 2 A)} 1+ w 2 ) 
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where c and d are constants which depend only on the dimension |H] . This is nothing but 
the universal 1/r decay found in the driven lattice gas and other related models [§J ll(Jj. 

It is now apparent that the most important quantity for determining the nonequilib- 
rium correction is the magnitude Q = c(±ei — > o) — c(o — > ±ei) of the quadrupole. By 
using (J3J) and the local detailed balance condition (|2|L we get 

Q = (1 - e ~ £ ) $(0, u; E) - (e E - 1) $(0, v; —E). (10) 

We of course have Q = for the exponential rule. For the heat bath rule, Q ~ vE 2 /2 in 
the lowest order. As for the Metropolis rule, we have Q ~ E 2 if \v | > and Q — f-E if 
|f | < again in the lowest order. 

Exact solution: To see that the power law decay is not an artifact of the perturbation, 
we show how one can get the exact solution of (JHJ). Our strategy is to define the effective 
charge q z which are nonvanishing only for z G {o} U U so that the exact solution ip z of 
(JHJ) also satisfies the Poisson equation J2 Z ' T ZjZ 'i/) z > — ~Qz- This is always possible since (JHJ) 
reduces to Y^ z < T Z)Z *if) z i = for any z G" {o} U U, and there are enough degrees of freedom 
in q z to fit the solution. Note that periodic boundary conditions and the (lattice version 
of) Gauss' law implies J2z Qz — 0. Then, by examining the symmetry of the model, we 
find that q z can be written with some Qi and Q 2 as q a = 2Qi + 2(d — l)Qi, q± ei = —Qi, 
and q~g = —Q2 for 5 G W . Therefore the asymptotic behavior (JUJ) is still valid if we replace 
Q with Qi - Q 2 . 

Note that ip z can be expressed as a linear combination ip z = J2z' Qz' G(z — z'), where 
G(z) is the solution of J2 Z ' T z>z 'G(z') = —5(z) where 5 Z = 5 Z)Q — L~ d . By substituting the 
Poisson equation J2 Z ' ^z,z'Vv = ~Qz into (JHJ), we get — q z + J2 Z ' D z , z >ip z > = —q z . This can 
be regarded as simultaneous linear equations for determining Qi and They can be 
solved, and one gets lengthy formulae (which we do not write down here) that involve 
G(z) as well as Q, s, s', s", t, and t' . In the lowest order, one of course finds Q\ ~ Q and 
Q2 ^ 0. 

Standard DLG: Finally we briefly discuss the case of the "standard" driven lattice gas 
jniUJIHl with hard core repulsion and nearest neighbor interaction. We restrict two-particle 
configuration (x,y) to those satisfying x 7^ y, and set V(x — y) — — J if \x — y\ = 1 and 
V(x — y) = otherwise. The analysis is essentially the same as the soft core case, and one 
again gets (J7J) as the condition for the absence of nonequilibrium corrections. We also get 
the long range correlation Q, where the "charge" Q should be "renormalized" jllj . But 
the basic behaviors are the same, and one has Q ~ —const. JE 2 in the heat bath rule, and 
Q ~ const. E 2 if \ J\ > \E\ and Q ~ — const. JE if \ J\ < \E\ in the Metropolis rule. 

Discussions: The universal l/r d long range correlation is often referred to as one of the 
essential features of nonequilibrium steady states in driven systems with a conservation 
law. The standard theoretical expression in terms of the fluctuating hydrodynamics jlj 
QUI may suggest that the long range correlation is a hydrodynamic or a many-body 
phenomenon. We have shown, however, that the driven lattice gas with only two particles 
may exhibit the same long range correlation, suggesting that it is essentially a two-body 
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phenomenon. It is quite likely that long range correlation in the two-particle model carries 
over to a many-body system. This is also suggested by a recent systematic expansion for 
the steady state of driven lattice gases [TT] . 

Perhaps more importantly, we have observed that the magnitude (or even the exis- 
tence/absence) of the long range correlation depends crucially on the rule of transition 
rates that one uses. It is vanishing i n a large class of models with the exponential 
condition (j7J), while it is fairly large (probably, too large) and may appear in the first order 
of E in models with the Metropolis rule. The disagreement is too sever to be regarded as 
a minor quantitative difference. We should probably start seriously examining which rule 
(if any) provides us with truly physical descriptions of nonequilibrium steady sates. We 
stress, however, that it is impossible to decide which rule is "realistic" by purely theoret- 
ical arguments, since we still do not know of any fundamental principles that determine 
nonequilibrium steady states. It would be desirable if one can find (just as in the case 
of equilibrium dynamics) some macroscopic constraints from a phenomenological point of 
view [HI] which enable one to determine microscopic dynamics |T1]. 

It is a pleasure to thank Joel Lebowitz, Raphael Lefevere, Elliott Lieb, Shin-ichi Sasa, 
Herbert Spohn for useful discussions. 
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